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Abstract 

We study the low-energy collective excitations and dynamical response functions of weakly cou- 
pled random antiferromagnetic spin-1/2 chains. The interchain coupling leads to Neel order at 
low temperatures. We use the real-space renormalization group technique to tackle the intrachain 
couplings and treat the interchain couplings within the Random Phase Approximation (RPA). We 
show that the system supports collective spin wave excitations, and calculate the spin wave veloc- 
ity and spectra weight within RPA. Comparisons will be made with inelastic neutron scattering 
experiments on quasi-one-dimensional disordered spin systems such as doped CuGe03. 



Antiferromagnetic (AF) quantum spin chains have been of interest to physicists since the 
early days of quantum mechanics^]. The one-dimensional nature of such systems allowed 
for tremendous theoretical progress both in clean systems by using exact solution and field 
theory mapping ]} Q] and disordered systems within renormalization group framework 0, U, 
S, 0, Q, 0). While such one-dimensional models have remarkably rich physics, in general 
they do not give a complete description of real systems. Real spin chain compounds, such 
as CuGeC>3 l| and KCUF3 lo| . always have some weak interchain couplings present, which 
can change the physics at lowest energy/temperature. For example, strictly one-dimensional 
models do not exhibit phase transitions into states with broken symmetry, while real spin 
chain systems often develop Neel order at very low temperatures due to the weak (3D) 
interchain couplings. It is thus important to study the effects of these interchain couplings 
to fully understand the low-energy/temperature physics of real spin chain compounds. 

In this paper we study the low-energy collective excitations and dynamical response 
functions of weakly coupled, disordered AF spin- 1/2 chains. Our work is motivated in 
part by the experimental studies on doped CuGeOs. In the absence of doping, it is a 
spin-Peierls system in which the spins dimerize and form a gapped, non-magnetic ground 
state. Upon doping, the system becomes disordered, and both dimerization and spin gap get 
suppressed. Amazingly, when doping reaches certain level the spins become Neel ordered 
at low temperature, which has been observed experimentally in Zn- and Si-doped CuGe03 

Since these experimental discoveries a number of 
S^JSWjWrtto* the static Nee! orde ri „ g in these HHH 

using mean-field theory. On the other hand the collective 
excitations and dynamical response functions, which have been studied experimentally using 
inelastic neutron scattering 20] , have not been studied theoretically thus far. The collective 
excitations and dynamical response functions are the subjects of the present work. We 
go beyond mean-field theory by allowing the Neel order parameter to fluctuate, and treat 
the interchain coupling using the random phase approximation (RPA), while tackle the 
intrachain coupling using the real space renormalization group (RSRG) method^, Q]. The 
RSRG technique has been proven to be powerful in obtaining ma gnet ic and thermodynamics 
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properties of random spin chains. Various numerical techniques 
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23| have also been 



deployed to study random chains and their results agree with those obtained by RSRG. 
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Our strategy here is similar to that of Schulz [3j|, who studied weakly coupled pure chains. 
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We find that despite the presence of disorder, the Neel state supports linearly dispersing 
spin waves, in agreement with experiments. We also obtain the spin wave velocity and the 
spectra weight of spin waves in the dynamical response function in terms of microscopic 
parameters of the system; this allows for detailed quantitative comparison between theory 
and inelastic neutron scattering experiments in the future. 

Consider weakly coupled spin- 1/2 antiferromagnetic chains with z nearest neighbor for 
each chain. The suitable Hamiltonian to describe this system is given by 

H = ^ Ji.riSifi • Si+irf + J± ^ ^ Sirf ■ Si rf+g, (1) 

where i is the site index along the chain, n is the chain index, and 5 is the index summed over 
the nearest neighbors. The intrachain couplings are drawn from a random distribution 
function P(Ji t H) (but with independent realizations for each chain), while the interchain 
coupling J± is taken to be constant. Both of the intrachain and interchain couplings are taken 
to be positive. Treating the interchain couplings in the mean-field appr oximat ion |36l. \!?\ , the 
presence of these couplings is effectively replaced by a staggered field which is responsible for 
long-range ordering at low temperature. The mean field approximation for the interchain 
couplings can be described as follows : for a given site i, the staggered field acting on this 
site is determined by the magnetization of neighboring sites siting on neighboring chains. 
In general the staggered field resulting from averaging the magnetization of the neighboring 
sites will be random. However in the limit of infinite coordination number z (say in the 
limit of large dimensionality), the number of neighboring sites contributing to the staggered 
field at site % becomes infinitely many and the fluctuations are suppressed; the staggered 
field becomes uniform in this limit because it becomes the average of the magnetization 
over infinitely many neighboring sites (see below). With this simplifying approximations at 
hand, the original problem is reduced to a random AF spin chain in the presence of uniform 
staggered field 

h 1d = J2 J * s i ■ s w - /^(-iysf. (2) 

i i 

The staggered field h is obtained through mean-field self-consistency condition 

h = zJ±m, (3) 

where m is the disorder-averaged staggered magnetization. The staggered magnetization 
at site i is rrii = (-iy(Si). The resulting effective ID problem can be solved using the 



real-space renormalization group (RSRG) method 



36). 



from which the phase diagrams of 



the systems have been obtained for various cases 

In the present work we go beyond the static mean field approximation and calculate the 
dynamical response functions by treating interchain couplings within the Random Phase 
Approximation (RPA), from which we also obtain the collective mode spectrum of the 
system. Within RPA the dynamic susceptibility is given byj^]: 

x a/3 

XRPA = " j - a /3 ( 4 ) 

where Xid * s the disorder-averaged single chain susceptibility matrix in the presence of the 
staggered field that satisfies Eq. (j2J). The expression given above is valid for transverse and 
longitudinal dynamic susceptibility. In this work we focus on the transverse response func- 
tion since the transverse part couples more directly to the collective excitation of the system 
than the longitudinal part does. Further motivation to study the transverse susceptibility 
is provided by recent experiments which focus on the transverse part; hence, a(3 = H — 
in Eq. (J1J. It is worth noting that although we concentrate on the transverse dynamic 
response in our calculation, the formalism developed here can be readily applied to obtain 
the longitudinal dynamic response and to study other random spin chains. 

Let us continue our discussion on the transverse dynamic response. For a specific disorder 
configuration, the chain susceptibility Xid represents the dynamical response of the chain at 
wave vector q' to an external perturbation at wave vector q; in general q 1 can be of any value 
due to the presence of disorder which breaks the translational symmetry. The symmetry is 
restored however once disorder averaging is performed (in fact the system is self-averaging), 
except for the doubling of unit cell by the staggered field. Thus a perturbation with wave 
vector q also induces response at another wave vector q' = q + tt, in addition to the usual 
response at q' = q. Hence the chain transverse susceptibility matrix is represented by a 2x2 
matrix : 

As a consequence of this, we can rewrite the RPA susceptibility, Eq. (jlj) as follows : 



xi D (q, 



— i / \ 1 / Xn(q,u) Xi2{q,u) , 

Xrpa(Q^) = . . . . ( 5 ) 
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where D(q,u) is the determinant of 1 — zJ±Xid (?> u )- Thus the calculation of Xrpa(Qj u ) 
reduces to the calculation of Xid(<1, cj )- 

To calculate Xid (?> s'> Ll 0> we use the RSRG method. In the RSRG scheme, one assumes 
the two spins that are coupled by the highest-energy bond form a ground state on this bond 
(singlet when there is no field); this bond is eliminated and new effective bonds between 
the remaining spins are generated perturbatively; this process is repeated until the ultimate 
low-energy limit is reached. Within this scheme the single chain response functions are 
approximated by the sum of those strongly coupled bonds that form during the RG process, 
and the coupling among different pairs of spins are neglected, as outlined in Ref. LiSj This 
approximation is asymptotically exact in the low-energy limit. The difference between the 
present case and that of Ref. I38II is that in addition to the AF bonds we also have the 
staggered field, which complicates the RG process. However it has been shown earlier [36J 
that in the limit of weak field (corresponding to weak interchain coupling), its effect on the 
RG flow is negligible and we thus do not consider it here. We thus start by considering a 
spin pair connected by a strong bond in the presence of a staggered field whose Hamiltonian 
is given by : 

H = ttS 1 -S 2 + h(S z 1 -S z 2 ), (6) 

where Q is the bond connecting the spin pair which we will identify as the cutoff of the 
system at a given stage of RG and h is the uniform staggered field as defined in Eq. (jHj); 
the eigen states of this Hamiltonian are : 

|o> = 1/^/1 + cg(c | ++> + !--» 
|i> = 1/1/1 + cf(c 1 | ++) + |--» 

|2> = | + +); |3> = |--> (7) 
with the corresponding energy eigenvalues : 

e = -n/4 - y/(n/2y + h*, 

E 1 = -fl/4 + v^/2) 2 + h 2 , 

E 2 = E 3 = fi/4, (8) 

where the coefficients C = 2h/Vt-^l + {2h/Q) 2 and d = 2h/Vt+^l + {2h/Q) 2 . Equipped 
with the spectrum of the spin pair, we proceed to calculate the dynamic structure factor for 



the pair. The T = spectral representation of the dynamic structure factor for a spin pair 
in the presence of uniform staggered field is : 

S&fabfc,") = ^<0|^Jm)(m|S-|0>%- AE) 
in 

= (c ° +e ;;g +e,m v - n /2 - vwn?) 

(9) 



where S^ 1 = S x + e iql S^ is the Fourier transform of the spin operator, I is the distance 

between two spins, |0) is the ground state and |m) are the excited states | + +) and | ) of 

such pair as written explicitly in Eq. ((7|). The system can only be excited to states different 
m z value because the operator connects states with different m z value; the transition 

can only happen from the ground state to the states | + +) and | ). 

To calculate the dynamics structure factor of the whole chain, we use theioint distribu- 
tion function of bond-length and strength, characterized in detail in Ref. y, and sum up 
contributions of all the strongly coupled bonds that are formed through the RSRG process. 



We follow the procedure outlined in Ref. |38| closely to obtain the dynamic structure factor 
of a single chain by summing the contribution from strongly coupled pairs; the dynamic 
structure factor for a single chain is given by : 

«!><*>) = n Tn J dld(P((,l;T n )S+ a ; r (q 1 ,q2,uj) 

(10) 

where P((, /; Tn) is the joint distribution of bond-length and strength. We follow the def- 
initions outlined in Ref. |4| to denote nr n as the fraction of spins left at energy scale Tn, 
( = ln(f2/ J) as the log energy scale, Tq = \og(fl /fl) as the log-cutoff of the energy scale, f2 
as the non-universal energy cutoff of the original Hamiltonian and Q as the energy cutoff of 
the renormalized problem. The transverse susceptibility for a single chain is then obtained 
by integrating the dynamic structure factor: 

where gi, q-i — q or q + 7r and S^(qi, q 2 , oj) is given by Eq. (JTHJ). 

We would like to study the collective excitations of the system, whose spectrum is given 
by the singularity of Xrpa^ EQ- ®> or the condition D(q,u) = 0. The rotational invariance 



of the system combined with the self-consistent condition Eq. (jSJ) requires D(q = 0,u = 
0) = 0, because the staggered field can be rotated without affecting the self-consistency. 
As a consequence the energy of the collective mode vanishes as the wave vector q goes to 
zero. Thus to obtain the collective mode dispersion it is natural to expand the quantities 
Xid (9>9> w )' Xiij(?> 1 + 7r >^)>XLD(<? + 7r >9>^)> and + 7r, g + 7r,a;) for small u and q: 

xlflfe?^) ^ xt D (0,0,0) + a g: qUJ 2 + b qA q 2 

XiD( ( i + 7r ^i + 7r ^) - xJflK^o) + %»,?+y 

- (12) 

where the expansion coefficients are given by : 

_ , du' S+n(0,0, a/) 

Or, 



' 7T d 



2k /"rib/ Co 1 a/ 2 + fr 2 

dw' 5+ D (0,7T,o/) 



7T d/ 2 



4 r dw' C 1 1 cj 2 + /i 2 

q ' q+n ~3J ~ 1 + C 2 ln(n /fi) ^ ^' 2 - /> 2 



where 



5+7(0,0,0;) 
5+ D (0,7r,a;) 
5+ D (7T,7r,w) 



(1 + Cp) 2 1 1 u 2 + h 2 
l + C 2 l v h?{VL /VL)uJUJ 2 -h 2 
l-C 2 1 lco 2 + h 2 
1 + C 2 l v \n 3 (n /Q)uuj 2 -h 2 
(1 - Cp) 2 1 1 to 2 + h 2 
l + C 2 ljn 3 (n /n)ujuj 2 -h 2 ' 



(14) 



where l v = a/(ln(fio/J)) is the microscopic length scale determined by the initial bond 
distribution; we use x to denote the variance of x. Using the condition that D(q,u) = 
det[l — zJ±Xii)(q, u))] = for uj = q = 0, we obtain quartic equation in ui 

Au* + Buj 2 + C = 0, (15) 



where the coefficients A, B, and C are : 



A = (zJ\ ) 2 a n „a 




B = -zJ±(l - zJ ± Xid (0,0,0))a g+7rj(?+7r 

+ {zJ ± ) 2 ({a q+W)g+w - a qA )b q , q q 2 - a 2 qq+7t ) 
C = (zJ ± (l-zJ ±X lD^AO))b q , q 

- (zJ ± )% q+n )q 2 - (zJ ± ) 2 bl q q\ (16) 

The solution to this quartic equation gives us the spin-wave dispersion of the system. To 
the leading order of the wave vector q we obtain a linear dispersing spin-wave us = v s q, with 

zJ ± m\n 3/2 (n /zJ ± m)l v . (17) 

yu 

In obtaining the result above for the spin-wave velocity, we have explicitly worked in the 
limit Jj_ — ► 0. 

We now calculate the dynamic structure factor within RPA, which can be accessed 
through Inelastic Neutron Scattering (INS) experiment. To obtain the dynamic structure 
factor near q = 0, we take the imaginary part of the upper left component of the RPA 
susceptibility matrix, Eq. (JHJ), i.e. Srpa = ^(xii(Qj a; ) / D(q, lj)) , 

where a small imaginary part, i5 has been introduced to shift the pole in the determinant 
D(q, uj) to slightly above the real axis. The pole in the determinant is realized when u = v s q, 
where v s is the spin- wave velocity given in Eq. (j!7j) . A straightforward calculation results 
in a simple form of the RPA dynamic structure factor near q = : 

Srpa = \^^5(u - v s q). (19) 

Following the same procedure, we also obtain the dynamic structure factor near q = it + 5q, 
Srpa = Ss(X22(q,w)/D(q,u)). The result is as follows : 

I zJ ± m 3 In 2 (VL /zJ ± m) 
Srpa = ^ 6{u, - v s (5q)), (20) 

where m is disorder-averaged staggered magnetization and Qq is non universal cutoff for 
the chain. Our calculation predicts that a sharp peak develops at the pole where to = v s q. 
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The intensity of the peak is proportional to the wave vector q(l/5q) near q = 0(tc). Within 
the Random Phase Approximation framework we only get a sharp peak of the intensity, 
as is shown by the delta function in the dynamic structure factor. The peak shows more 
pronounced contribution from the dynamic structure factor near q = tc because the long- 
range staggered configuration is realized near this wave vector. We comment here that 
while our expression for v s (Eq. ()17|) involves parameters of the random distribution (Qq 
and l v ) that cannot be measured directly, our results on the spectral weight can be compared 
directly with the intensity of inelastic neutron scattering experiment, once v s is determined 
from the measurement; this is because Eq. ()19j) involve v s and other measurable quantities 
only, and the same combination of Qq and and other measurable quantities appears in Eq. 
ffT7j) and (J2TIJ) (l v is of order one lattice spacing for generic distributions). Thus our results 
allow for detailed quantitative comparison with future experiments. 

A few years earlier, Martin et al. studied the excitation spectrum of doped CuGe0 3 us- 
ing inelastic neutron scattering^, and found sharp propagating spin wave excitations when 
the system is Neel ordered, despite the fact that the Neel phase is stabilized by disorder. 
They found the spin wave spectrum to be linear. Our results agree with these experimental 
findings, and it is clear that such propagating excitations must be collective modes sta- 
bilized by interchain couplings, as single random chains do not support such propagating 
modes|38[. In the present work we assume there is no dimerization, while in doped CuGeOs 
dimerization survive and coexist with Neel order when doping level is sufficiently low. It 
is straightforward to generalize the present approach to the case with dimerization j^, Si. 
as well as finite temperature and chains with other spin sizes. Recently Masuda et al. [39] 
studied the dynamic spin-spin correlation of a new compound BaCu 2 (Sii_ x Ge a; )207 using 
inelastic neutron scattering. This system can be described very well by antiferromagnetic 
spin-1/2 chains with random exchange due to the random distribution of Si and Ge atoms. 
The experimental data on the d yna mic structure factor on this compound fit the universal 
scaling form predicted in Ref. |3j| very well. Our theoretical work proposed here could 
be of relevance to this new experimental realization of random exchange antiferromagnetic 
spin-1/2 chain; in particular it would be interesting to study the collective excitation on this 
compound in the ordered phase and compare it with the results obtained here. We hope the 
present work will motivate future experiments that will study the spectral weights of the 
spin waves in detail, and test the predictions made here on them. 
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